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We consider two recently proposed measures of non-Markovianity applied to a particular quantum process 
describing the dynamics of a driven qubit in a structured reservoir. The motivation of this study is twofold: on 
one hand, we study the differences and analogies of the non-Markovianity measures and on the other hand, we 
investigate the effect of the driving force on the dissipative dynamics of the qubit. In particular we ask if the 
drive introduces new channels for energy and/or information transfer between the system and the environment, 
or amplifies existing ones. We show under which conditions the presence of the drive slows down the inevitable 
loss of quantum properties of the qubit. 
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I. INTRODUCTION 

The inevitable interaction of a quantum system with its en- 
vironment leads to dissipation of energy and loss of quantum 
coherence 01 lit- Unfortunately these phenomena can funda- 
mentally limit the potential of quantum technologies, whose 
power is based on quantum effects [3]. In suitable surround- 
ings, however, the quantum system may temporarily regain 
some of the previously lost energy and/or information due to 
non-Markovian effects in the system dynamics. This is one 
of the reasons why the study of non-Markovian quantum dy- 
namics has received an increasing amount of interest in the 
last several years H. It is surprising, then, that the definition 
of non-Markovianity still remains elusive and, in some sense, 
controversial. 

Markovianity is well defined for classical stochastic pro- 
cesses: a stochastic process has the Markov property if the 
probability distribution for the future states of the process de- 
pends only on the present state. Loosely speaking, depen- 
dence on past states should then be a characteristic trait of 
non-Markovian processes. The mathematical definition of 
non-Markovianity does not, however, translate easily into the 
widely used language of density matrices and master equa- 
tions in quantum physics. In the case of an open quantum sys- 
tem evolution described by a Lindblad master equation 
one can find a stochastic descriptions for the system dynam- 
ics, e.g., with the Monte Carlo Wave Function method, where 
future dynamics of individual quantum trajectories only de- 
pend on their present states BHD. In this spirit some physi- 
cists have attributed non-Markovianity to the breakdown of 
the semi-group property ||9[] or a generalization thereof, the 
divisibility property [10]. However, finding a stochastic de- 
scription that corresponds to more general master equations is 
not, in general, a straightforward task 111 ill . 
With the aim of extending the concept of non-Markovianity to 
more general quantum dynamics, efforts have been made to 
clarify the very definition of non-Markovianity in the context 
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of open quantum systems and to quantify the amount of non- 
Markovianity in a given quantum process. One viewpoint as- 
sociates non-Markovianity to dynamical dependence on past 
evolution. This may mean either dependence of the evolution 
on all past states of the system, as in the memory kernel ap- 
proach Il2ill3ll . or a dependence of the asymptotic state of the 
system on its initial state fl4ll . Yet another definition equates 
non-Markovianity with a partial recovery of previously lost 
information IU5l - [l7ll . It is worth stressing that the definitions 
of non-Markovianity are generally not the same, nor are they 
mutually exclusive, and that they may agree perfectly about 
the non-Markovian character for some models and disagree 
about others lfl8ll . 

In this paper we consider two recently proposed measures 
of non-Markovianity introduced by Rivas, Huelga and Ple- 
nio (RHP) OH and Breuer, Laine and Piilo (BLP) 0] and 
apply them to a quantum process describing the short and in- 
termediate time-scale dynamics of a driven qubit in a struc- 
tured reservoir niHU. The comparison allows, on one hand, 
to elucidate the differences and analogies between the non- 
Markovianity measures in a physically interesting model, and 
on the other hand to gain new understanding on the effect 
of the driving force on dissipative dynamics. For the sake 
of concreteness in the rest of the paper we will refer to the 
case of a laser-driven atom. All the conclusions, however, are 
valid for any driven two-state system. Due to the presence 
of the drive, the commonly used secular approximation does 
not always hold, especially in the short non-Markovian time 
scales. The importance of the nonsecular terms has been re- 
cently demonstrated in the context of superconducting circuits 
inRef. JUIH. 

The article is organized as follows. In Sec. [II]we introduce 
our model, namely the driven qubit, and the master equation 
describing its dynamics in a structured environment. In Sec- 
tions [III] and [IV] we introduce and study the RHP and BLP 
measures of non-Markovianity using this model. In Sec. [V] 
we compare the non-Markovian properties of a driven qubit 
with those of a non-driven qubit and finally, in Sec. [VI] we 
summarize the results and present our conclusions. 
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II. THE DRIVEN QUBIT 

We consider a qubit with energy separation uia (fi = 1), 
driven by a laser of frequency lu^ almost resonant with the 
transition frequency of the qubit: |A| = \uia — u>l\ <C u>a- 
The qubit is embedded in a zero-T environment. The strengths 
of the interaction between the qubit and laser and the qubit 
and environment are quantified by the Rabi frequency 57 and 
the coupling constant a, respectively. We assume that ft <C 
uja, a condition that is typically satisfied in quantum optical 
situations. When a is the smallest relevant frequency, that 
is, the qubit couples weakly to the reservoir, the dynamics of 
the qubit is given by the following local-in-time completely 
positive master equation, derived by one of us in Ref. 11911 : 



dp(t) 
dt 



-i[H, pit)} + T> s (t)p(t) + V NS {t)p(t), (1) 



where the unitary evolution of the qubit is generated by the 
system Hamiltonian 
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Va 2 + o 2 , 



(2) 



and the dissipative terms have been divided into a secular con- 
tribution T>s and a nonsecular contribution T>ns- Introducing 
the short-hand notation 

C[A(t)Mt)} = 7(*) (A(t) P {t)AHt) - l -{AHt)A{t), P {t)} 

(3) 

to describe a Lindblad-type decay channel with a jump opera- 
tor A(t) and a corresponding decay rate -f(t), we can express 
the dissipative terms as: 



V s {t)p(t) = ClC[a_ rl+ {t)] + C 2 _C[a +r i^t)] 
+ Cgr[o-„7o(t)], 



(4) 



V NS (t)p{t) = T-(t){C-C [a +P (t)<j z - a z a+p(t)] 

+ C+C- [a+p(t)a+ - a+a+p(t)] } 

+ r + (i){c+C [o-p(t)o z - a z a_p{t)\ 

+ C+C- [<r-p(t)tr- - <T-a-p{t)) } 

+ r (t){c.C [a z p(t)a. - a-a z p(t)] 

+ C+C [a z p(t)a+ - a+a z p{t)} } + H.c.(5) 

where C± = (A ± u>)/2u>, Cq = ft/2u> and H.c. denotes 
Hermitian conjugation. The master equation is given in the 
dressed state basis |^±) = ± J C± | e) + J C T \g), where |e) 
and \g) are the excited and ground states of the qubit. The 
decay rates in T>n$ and T^s are connected via the relation 



r«(t) 



7«(*) 



-iX^t) 



(6) 



where 7^, € R and £ G {—,0, +}. For the sake of con- 
creteness we focus on a qubit embedded in a reservoir with a 



Lorentzian spectral density J(ui) — a/(2ir) A 2 /[A 2 + (ui — 
Wo) 2 ], where uio is the center frequency and A the width of the 
Lorentzian. For this case the decay rates in Eqs. (@])-(|5]l take 
the form JH 



7«(r) 



2(1 + g a ) ( 1 ~ 6 Tcos 1t T + e T ^sin^T) 



A C ( T ) = TT-^ (-96 + e cos 9( T + e " T sin . 

1 + 

(7) 

where T = Xt and = s — £p, again with £ = {—,0, +}. 
Different dynamical regimes are defined in terms of the pa- 
rameters 



X : 



T C UJ 
T S X 



(8) 



The parameter s quantifies the detuning between the qubit 
and the center frequency of the Lorentzian. The parameter 
p expresses the relationship between the reservoir correla- 
tion time-scale tc = A -1 and the typical system time-scale 
ts = ui^ 1 and it determines the border between secular and 
nonsecular regimes. Note that here we refine the typical text- 
book notion of the secular approximation, which requires that 
the typical time-scale of the system is negligible in compar- 
ison with the relaxation time-scale and results in a coarse- 
graining of the relaxation dynamics. Since in this work we 
are interested in the short time-scale dynamics we assume a 
stronger condition, namely that the typical time-scale of the 
driven qubit is much smaller than the other time-scales and, in 
particular, the reservoir correlation time-scale tc- We call the 
limit ts <C tc the strong secular limit. 



A. Strong secular limit 

When p ^> 1 the secular approximation holds and the non- 
secular dissipation term 2?jvs can be neglected from Eq. (fl}. 
Consequently, in the language of Markovian |7j, l8|] and non- 
Markovian 1I23II24I1 quantum jumps, the master equation com- 
prises of three Lindblad-like terms describing phase flips and 
jumps between the eigenstates of the driven qubit, with the di- 
rection of the jump (regular/reversed) depending on the sign 
of the corresponding decay rate (positive/negative). We have 
shown in Ref. [20] that, in the secular regime, the decay 
rates j±(t) always oscillate taking temporarily negative val- 
ues, whereas 70 (t) is positive for small values of s and has 
periods of negativity when s > 3.6. 



B. Nonsecular limit 

In the opposite regime, defined by p -C 1, we have to con- 
sider the full master equation including the nonsecular terms 
to obtain a proper description of the short time-scale dynam- 
ics. However, in this limit the decay rates take a simpler form. 
More precisely, 7±(t) ~ 7o(*) = l(t) and X±(t) Xo(t) = 
X(t) and again j(t) is positive for small values of s and takes 
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temporarily negative values when s > 3.6, while A(i) < for 
all times t. With this approximation the master equation can 
be cast into a remarkably simple form with a single Lindblad- 
type decay channel: 



dp(t) 
dt 



= -i[H + H',p(t)]+£[A ) j(t)], 



(9) 



where the jump operator is 

A = C-<T + +C+<T-+C <T Z , (10) 

and the additional term in the coherent evolution is 

H' = \(t)C (C + -C_)(<T_+a + ). (11) 

III. SEMI-GROUP AND DIVISIBILITY 

Historically non-Markovianity has been closely associated 
to deviations from the Lindblad master equation. The com- 
pletely positive trace preserving (CPTP) map $ : p(0) >-> 
p(t) = Q(t)p(0) corresponding to a Lindblad master equation 
always satisfies the semi-group property $(i+s) = <l>(i)<i>(s). 
For example Wolf, Eisert, Cubitt and Cirac propose to call a 
map Markovian if it is a CPTP map satisfying the semi-group 
property |0] . The master equation (Q]l for the driven qubit in 
any structured reservoir is not in the Lindblad form due to the 
time dependency of the decay rates and therefore the dynami- 
cal map is not an element of a single parameter semigroup and 
Markovian in this sense. 

The measure for non-Markovianity proposed by Rivas, 
Huelga and Plenio is based on a more general property of 
maps than the semi-group property, namely divisibility 11X011 . 
A CPTP map + r, 0) is divisible if it can be written as a 
decomposition of two CPTP maps $(£ + r, t) and 0): 



$(i + r,0) = *(* + r,i)$(t,0). 



(12) 



When the dynamical map is homogenous in time, &(t + 
r,t) = $(t, 0) = $(r), Eq. (fLZb reduces to the semi-group 
property. 

RHP define a map to be Markovian exactly when it is divisi- 
ble. This amounts to requiring that <£>(£ + t, t) is completely 
positive. It is shown in Ref. fllOll that the quantity 



g(t) = lim 



|*(t + e,t)®I|0)M||-l 



> 



(13) 



is strictly positive if and only if $ is indivisible and hence g(t) 
identifies non-Markovian processes, according to their defini- 



tion. Here \<t>) — Sj=o K)l*)> wnere d is the dimension 

(|11> + |00»/V2 



of the quantum system; in our case 
is the maximally entangled Bell state. The quantity g(t) gives 
rise to a measure for non-Markovianity defined as 



dsg(s), 



(14) 



with the property that Nrhp G [0, 1]. The lower and 
upper limits correspond to Markovian and maximally non- 
Markovian maps, respectively. 



For a local-in-time master equation dp(t)/dt — C(t)p(t) the 
dynamical map, in the limit e — > 0, is formally given by 
3>(t + e, t) — exp[C(t)e\. Consequently, RHP propose that, in 
this case, one may check the non-Markovianity of the dynam- 
ical process by means of a simplified quantity 



9 (*) = I™ 



||{l + e[£(i)®I]}|0>M|| -1 



A. Secular regime 



(15) 



A direct calculation of Eq. (fl~5T > in the case of a driven qubit 
in the strong secular regime gives 



C 2 + P[l + {t)] 



ClP[ 1 -(i)]+2CtP[ l0 (t)\ 



,(16) 



where we define an auxiliary function P: P(x) = when 
x > and P(x) = —x when x < 0. In general, when- 
ever any of the decay rates is negative, divisibility is lost. In 
the Lorentzian case this is, indeed, the situation in the secu- 
lar regime where j± (t) < always for some periods of time. 
Thus we prove that according to RHP the dynamics of the 
driven qubit in the secular regime and in the Lorentzian case, 
is always non-Markovian. 



B. Nonsecular regime 



In the nonsecular regime we find 



g(t) = 



[C\ + CI+ 2Cb)P[ 7 (t)] 



(17) 



that is, the dynamics is indivisible and hence non-Markovian 
according to RHP whenever the decay rate j(t) takes tem- 
porarily negative values. Again the result is valid for any 
reservoir, in the weak coupling limit. Specifically, for a 
laser-driven qubit in a Lorentzian reservoir divisibility is bro- 
ken whenever s = (cjq — > 3.6, that is, the laser 
and the atom are sufficiently detuned from the center of the 
Lorentzian, as discussed in Sec. 



IV. INFORMATION BACKFLOW 

As a second measure of non-Markovianity we consider 
the one proposed by Breuer, Laine and Piilo in Ref. lfl5ll . 
The BLP measure aims at identifying non-Markovian dynam- 
ics with certain dynamical physical features of the system- 
reservoir interaction, rather than with the mathematical prop- 
erties of the dynamical map of the open system. In particular, 
BLP define non-Markovian dynamics as a time evolution for 
the open system characterized by a temporary flow of infor- 
mation from the environment back into the system. This re- 
gain of information manifests itself as an increase in the dis- 
tinguishability of pairs of evolving quantum states. Hence, the 
dynamical map $ is non-Markovian according to BLP if there 
exists a pair of initial states p\ (0) and pi (0) such that for some 
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time t the distinguishability of pi(t) and p 2 (t) increases, that 
is, 

a( Pl ,p 2 ,t) = -D[p 1 (t),p 2 {t)] > 0, (18) 

where D(pi,p2) — \\\pi ~ P2I is the distinguishability of 
pi and P2 with \\A\\ = V AA^ being the trace distance, and 
Pi/ 2 (t) — &(t, 0)pi/ 2 (0). When the distinguishability of two 
states increases information flows from the environment back 
to the system and vice versa. With this definition one can 
quantify the amount of non-Markovianity in the quantum pro- 
cess as follows 

Nblp{$) = max / dt a(p u p 2 ,t), (19) 

Pi,a(0) J a>0 

that is, Nblp gives the maximum amount of information that 
can flow back to the system for the given process. 
All divisible maps continuously reduce the distinguishability 
of quantum states, that is, for a divisible $ we get Nblp = 0. 
Hence the connection between the BLP and RHP measures is 
the following: if a map is Markovian according to RHP, it is 
Markovian also according to BLT. The converse is in general 
not true. 



A. Secular regime 



hold when the driven qubit is embedded in some other struc- 
tured reservoir. Indeed, negativity of one of the decay rates 
does not guarantee that a(t) > 0, although it immediately 
breaks the divisibility. It is therefore natural to ask if we 
can find some general parameter regions when the dynamics 
of the driven qubit is Markovian according to BLP, but non- 
Markovian according to RHP. It turns out that the question is 
non-trivial from many points of view. For one, showing that 
the reduced dynamics of an open system is non-Markovian 
according to BLP is simple in comparison to showing that 
the dynamics is Markovian, since the former case amounts 
to finding one pair of initial states whose distinguishability in- 
creases for some time interval, while the latter case requires 
that one probes all possible pairs of states and shows that dis- 
tinguishability can only decrease. Such optimization proce- 
dure is possible in practice only numerically, making it hard 
to isolate parameter regimes for different dynamical behav- 
ior. Moreover, the BLP non-Markovianity quantifier depends 
strongly on the spectral density of the environment via linear 
combinations of the three decay rates and their integrals, so 
identifying the parameters for which the two measures deviate 
may not even be possible unless one has a specific environ- 
mental spectral density in mind. Moreover, even specifying 
the spectral density, one needs an extensive numerical study 
involving several variables. So far we have been unable to 
find instances when the two measures disagree and the ques- 
tion remains open. 



A straightforward application of Eq. ( fT8l to the solution of 
the master equation (Q~|i in the secular regime gives 

e -mt)A'(t)(5x 2 + Sy 2 ) + e~ 2T ^T'{t)Sz 2 

c(t) = — , 

2x / e - 2A W(Sx 2 + Sy 2 ) + e- 2T ^Sz 2 

(20) 

where Sx — x\{Q) — #2(0), 5y — yi(0) — 2/2(0), and 
8z = z\ (0) — z 2 (0) are the differences in the x-, y- and z- 
components of the Bloch vector representations of p\ (0) and 
p 2 {G), respectively, and 

T{t) = lJo dS [ C + 7+(s) + °- 7 - (s) + 4C o7o(*)] , 

A(t) = f ds [C 2 +1+ (s) + C 2 _j-(s)] . (21) 
Jo 

The expression in Eq. d20l > is valid for a generic structured 
reservoir and the non-Markovian properties of the driven qubit 
strongly depend on the choice of the environment. We now 
focus on the Lorentzian case. Let us choose a pair of initial 
states such that 5z = 0. Then the sign of a only depends on 
A'(t) = C 2 + ~f + {t) + C^7_(*). Recall from Sec. ED that in 
the secular regime j± (t) have negative periods for all possi- 
ble parameter values. Then a(t) > always for some periods 
of time. Therefore, in the secular regime the process describ- 
ing the dynamic of a driven qubit in a Lorentzian reservoir is 
always non-Markovian according to BLP. 
It is worth stressing that although the BLP and RHP mea- 
sures both agree that the dynamics of the driven qubit in a 
Lorentzian reservoir is non-Markovian, the result need not 



B. Nonsecular regime 

In the nonsecular regime the master equation ([TJ is not, in 
general, solvable analytically and we can evaluate the quantity 
(7 only numerically. We can, however, consider the special 
case of resonance between the atom and the laser. In this case 
the master equation is solvable and we find 

l{t)e- 2 fo [( fa 2 _ g z 2) + 2Sy2] 

°W = y 

V2y(Sz 2 + Sx 2 ) + e- 2 Jo ds ~t( s ) [(Sx 2 - Sz 2 ) + 2Sy 2 } 

(22) 

Thus in the resonant nonsecular case the process is non- 
Markovian according to BLP if and only if 7(4) < for some 
time t. Numerical studies of the off-resonant case A 7^ 
indicate that this result holds more generally, that is, in the 
nonsecular regime A/"sLp(i) = if and only if 7(4) < 0. 



V. LASER INDUCED NON-MARKOVIANITY 

In this section we compare the non-Markovian character 
of the driven dissipative qubit and the analytically solvable 
model of an unperturbed dissipative qubit when both are em- 
bedded in a Lorentzian reservoir Mil 12511 . The master equation 
for the qubit in a Lorentzian reservoir, in the special case when 
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the qubit is in resonance with the Lorentzian, uja — wq, is 



dp(t) 
dt 



l(t) 



(23) 



where, in contrast to the master equation of Eq. ([]]), the opera- 
tors are in the bare state basis {|e). \g)}. The time-dependent 
decay rate is 



lit) 



2a\smh(dt/2) 



dcosh(dt/2) + Xsinh(dt/2)' 



(24) 



where d — \/A 2 — 2a\. Applying this model to Eqs. ([T8T i and 
( fT5T > and introducing quantity T(t) = J* dsj(s)we get 

j(t) [2e-^6z 2 + e-^(Sx 2 +5y 2 )] 
a \t) = . / „,„ . „ / ±\ / f „ s ( 25 ) 



4 v / e -2r(t) ( 5 z 2 + e -r(t)( 6x 2 + Sy 2) 



g(t) = -P[i(t)]- 



(26) 



Equations (125ll-(l26li show that, in absence of the driving laser, 
the dynamics is non-Markovian for both the RHP and the BLP 
measures if and only if ^(t) < for some time t 6 R + . This 
happens exactly when A < 2a, that is, when the qubit couples 
strongly to the environment. 

Consider now a situation where the qubit couples weakly to 
the Lorentzian reservoir and the transition frequency of the 
qubit is resonant with the center of the Lorentzian. The anal- 
ysis above shows that, without the laser drive, the dynamics 
of the qubit is Markovian according to both RHP and BLP. If, 
instead, we drive the qubit with a laser, the results obtained in 
Sees. II VI and Hill show that we can induce non-Markovianity 
in the system dynamics when the laser parameters are suitably 
chosen. More precisely, the laser should couple to the qubit 
in such a way that condition p ^> 1 holds, that is, we are in 
the secular regime. The secular regime can be achieved with 
a large enough Rabi frequency £1 and/or when the laser -and 
consequently the atom- are sufficiently detuned from the cen- 
ter of the Lorentzian. When this happens the laser interferes 
with the dynamics of the dissipative qubit so much that it can 
temporarily reverse the flow of information from the qubit to 
the environment. 



using two different non-Markovianity measures, and we have 
analyzed the effect of the driving laser on the non-Markovian 
properties of the dissipative qubit. 

In the secular regime both measures confirm that the dynam- 
ics of the driven qubit is always non-Markovian when it is 
embedded in a Lorentzian reservoir. However, in general this 
need not be the case and deviations between the two measures 
could be discovered when the driven qubit is embedded in an- 
other reservoir. So far we have been unable to find a spectral 
density for which the two measures disagree. This might be an 
indication that there are no realistic physical situation where 
this happens, even if the mathematical definitions of the quan- 
tifiers suggest this. It remains an interesting open question 
whether the two non-Markovianity measures always agree in 
the case of a driven qubit embedded in a structured reservoir; 
indeed deviations between the two measures are neither char- 
acterized nor well understood to date. 

Instead in the nonsecular regime the non-Markovianity mea- 
sures agree perfectly. In this regime the appearance of non- 
Markovian features is strongly dependent on the way the laser 
is coupled to the qubit. Neither the presence of the struc- 
tured reservoir nor the presence of the laser-drive is sufficient 
to guarantee non-Markovianity in the system dynamics in the 
nonsecular regime. In Sec. [V] we investigated the origins of 
these non-Markovian effects by comparing driven and non- 
driven qubits in a Lorentzian reservoir; it was shown that non- 
Markovian effects are not possible without a strong laser-drive 
and/or detuning between the driving laser and the center fre- 
quency of the Lorentzian distribution. This discovery has a 
clear physical interpretation; only a strong enough laser drive 
can induce a feedback of information from the environment 
into the system. In fact, non-Markovianity in the nonsecular 
regime occurs in the far off -resonant case. Far off -resonance 
dynamics is usually associated to the presence of virtual pro- 
cesses which may be at the origin, in this case, of the flow- 
back of information. 
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